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Simulation  of  a  35  GHz  Klystron 
Using  a  Magic2D  Template 


The  following  document  describes  the  steps  involved  in  tuning  a  6  cavity  35  GHz  klystron.  The 
initial  design  for  this  template  was  developed  at  SLAC  (Caryotakis,  Scheitrum,  Sphren,  and  Song). 
Modifications  have  been  made  at  MRC.  The  bulk  of  the  klystron  template  has  been  buried  in  an 
internal  MAGIC  file.  The  remainder  of  the  template  consists  of  the  controls  (type  of  simulation)  and 
the  tuning  parameters.  It  is  recommended  to  the  user  that  the  control  file  be  used  to  annotate  the 
results  of  each  simulation.  Since  the  tuning  procedure  is  iterative,  it  will  be  instructive  as  well  as 
essential  to  record  the  sequence  of  changes  within  the  control  file,  (or  within  a  sequence  of  modified 
files.  Each  set  of  changes  potentially  impacts  the  performance  of  each  of  the  following  simulations, 
so  good  record  keeping  is  critical. 

The  magic  input  file  for  the  klystron  template  is  “F3506.m2d”.  This  file  includes  the  user  modifiable 
input  and  tuning  parameters.  Once  these  are  set  the  file  calls  the  internal  coded  klystron  template 
embedded  in  MAGIC2D.  This  document  indicates  the  procedure  and  order  in  which  the  various 
tuning  parameters  should  be  used. 

At  the  beginning  of  the  file  are  the  current  design  parameters  for  the  Frequencies,  Qe,  Qo,  and  R/Q 
for  the  input  cavity,  the  intermediate  cavities,  and  the  output  cavity. 

The  file  is  setup  so  that  you  serially  run  a  series  of  Steps  beginning  with  Cavityl  (the  input  cavity) 
thru  Cavity6  (the  output  cavity. 

The  cavity  tuning  procedure  for  Cavities  1  through  5  is  the  following. 

1.  Run  one  of  the  Steps,  Step  1  thru  5. 

2.  Obtain  the  frequency. 

3.  Adjust  the  tuners  to  increase/decrease  the  frequency. 

a.  The  tuner  parameters  are  tune_left  and  tune jight. 

b.  If  frequency  is  high,  decrease  values  of  tune_left  and  tune_right.  Rerun. 

c.  If  frequency  is  low  increase  values  of  tune_left  and  tunejright.  Rerun. 

d.  If  frequency  is  good,  then  frequency  tuning  is  complete. 

4.  Check  Qeffective  and  R/Q.  Remember  that  you  are  measuring  and  tuning  to  Qeffective,  not  to  Qe 
orQ0. 

a.  If  the  measured  Q  is  greater/less  than  the  desired  value  of  Q  than  adjust  the  loading 
conductance.  Adjust  the  conductance,  sigma,  by  the  following  factor.  Sigma 
(desired)  =  Sigma  (used)  *  Q  (measured)  /  Q  (desired). 

b.  Rerun  and  verify  Q. 
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Selecting  Stepl  allows  you  to  tune  Cavityl.  The  section  of  the  file  for  this  Step  looks  like  the 
following: 


!****************************************************************************** 
!==>  Stepl,  Cavityl,  Can  Use  either  Magic2d__Sng  or  Magic2d_Dbl 
!  Cold  test  and  tune  cavity  1. 

I****************************************************************************** 

Comment  "  ===  Tuning  for  Input  Cavity  ==============:================== ======="  ; 

IF  (ISTEP.eq. 1)  THEN  ; 

I_Run__Type  =1;  !  0  =  Hot  test 

Cold_Test_Cavity  =  1  ; 

ENDIF  ; 

!  Cavity  1,  tuning  parameters. 

Parameter_GHZ  Desired_fl  =  fl;  !  Desired  Frequency  for  Cavity  1 

Tune_Left_l=  13  ;  !  Use  tuners  to  adjust  the  frequency, 

Tune_Right_l=  12;  I  Change  by  about  +/-1,  to  give  about 

+/-  6.5  MHz  frequency  change. 

PARAMETER  QT1  =  1/ (l/QEl  +l/Qol  )  ;  !  Desired  Q  effective  (of  Qo  and  Qe) 

Parameter  Sigma__l  =  0.32114E-01;  !  Change  this  value  to  adjust  the  Q. 

!  Record  results  here. 

ParameterJ3HZ  CF_1  =  34.982  GHZ  ;  !  Frequency 

Parameter  CQ_1  =  298  ;  !  Q  effective  (of  Qo  and  Qe) 

Parameter  CROQ  1  =  127  !  R  over  Q  as  measured  on  Axis. 


When  you  run  Stepl  you  should  see  a  figure  like  the  following.  This  illustrates  the  Ez  field  in  the 
cavity  and  the  beam  tunnel. 
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The  figure  of  interest  will  give  you  the  Q  of  the  cavity.  It  will  look  like  the  following. 

(Please  note,  that  are  multiple  measurements  of  Q,  but  they  should  all  give  the  same  value  within  a 
percent  or  so. 
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Using  the  mouse,  you  can  the  precise  value  of  Q  as  298.  Record  this  value  in  the  control  file  as 
CQ_1  =  298.  (See  above  listing.)  We  measure  Q  by  directly  looking  the  energy  in  the  cavity  and  the 
rate  at  which  it  is  decreasing  in  time.  Thus  we  use  Q  =  -2nf  W(t)/(dW(t)/dt).  In  the  above 
measurement  of  Q,  we  start  with  total  energy  in  the  simulation,  evaluate  the  derivative,  and  we  know 
the  frequency. 

The  next  figure  of  interest,  will  look  like  the  following.  Use  the  mouse,  click  on  the  curve  in  the  flat 
region.  Record  the  value  of  the  frequency,  cf_i  =  34.982  ghz.  This  is  the  resonance  value  of  the 
cavity. 


Signal  analysis  Frequency  for  E.DL  at  AXI  S_CAV_1 


The  next  figure  of  interest  will  look  like  the  following  this  a  measurement  of  the  R/Q  on  axis. 
(Incidentally,  we  also  have  such  a  measurement  at  the  cavity  gap  radius.)  Record  the  value  of  R/Q, 
as  croq_i  =  127.  To  obtain  R/Q  we  use  MAGIC  to  measure  the 

<VaxjS2>(t),  the  integral  of  the  Ez(z,r=0),  averaged  over  1  RF  period,  and  to  measure  <We>(t),  the 
electric  field  energy  in  the  cavity,  averaged  over  1  RF  period.  Then  the  ratio  of  these  measurements 
gives  us  a  measure  of  R/Q,  as  follows. 

<R/Q>(t)  =  <Vaxis2>(t),  /  (2nf  <WE>(t)). 

Initially  we  are  driving  the  cavity,  so  this  value  must  be  obtained  after  the  drive  is  terminated  the 
cavity  has  relaxed  to  its  natural  resonant  frequency. 
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TRANSFORM:  R  S  HUNT_T  C_Q^_AX  IS_1 


The  same  procedures  hold  for  Cavities  1  through  5.  These  are  the  single  gap  cavities. 


Run  Steps:  1  through  5  and  record  results  and  parameter  updates  in  control  file. 


Eample  Step  6:  Tuning  the  Frequency  of  the  Output  Cavity 

In  Step=6,  we  tune  the  Frequency  of  the  output  cavity.  Unlike  cavities  1  through  5,  we  use  the 
eigenmode  solver  to  get  a  good  estimate  of  the  frequency.  This  cavity,  which  has  two  gaps,  can 
oscillate  at  two  frequencies,  the  0-mode  (tune  to  35.00  GHz)  and  the  7t-mode  (which  is  about  0.5 
GHz  below  the  0-mode  frequency).  The  proximity  of  these  two  modes,  make  it  difficult  to  tune  by 
exciting  the  cavity  with  a  drive  signal,  since  any  drive  envelope  will  have  some  contamination  of  the 
undesired  component.  It  is  preferable  to  use  the  MAGIC2D_DBL  in  obtaining  this  solution,  as  we 
can  make  direct  measurements  of  R/Q  from  the  eigenmode  solution,  provided  the  fields  have  only 
weak  contamination. 

J  ****************************************************************************** 

!  s==>  Step6,  Output  Cavity,  Use  Magic2d_Dbl 
1  Cold  test  and  tune  Output  cavity. 

I  This  output  cavity  is  set  for  0-mode,  not  the  pi-mode. 

I  Output  Cavities,  USe  double  precision  Eigenmode  to  tune  these  cavities. 

1  Get  R/Q  FROM  EIGENMODE,  Integrate  |Ez|  along  axis. 

1  Then  R/Q  =  ( |e| .dl) A2/4pi/f /ljoule_from  eigen  soln. 
!****************************************************************************** 


IF  (ISTEP.eq. 6)  THEN  ; 

I_Run_Type  =  2  ; 

Cold_Test_Cavity  =  OutputCAvity  ; 
I__Cold_Type  =  2  ; 

ENDIF  ; 

Parameter_GHZ  Eigenmode_fOut_0  =  Fout_0  ; 
Tune_Lef  t__6=  8; 

Tune_Right_6  =  8; 

Tune  Left  7*=  8; 


t 


Tune_Right_7=  7; 

1  Record  Results. 

Parameter_GHZ  CF_7_0  =  34.998  GHZ;  !  adjusted  after  eigenmode. 

Notice  that  this  cavity  has  4  tuning  elements.  These  correspond  to  shims  in  each  side  of  the  dual 
cavity  output  structure.  By  adjusting  the  number  of  cells  in  each  shim  element,  the  frequency  can  be 
fine  tuned  to  within  6MHz. 


The  first  figure  of  interest  in  running  Example  Step  6  is  the  contour  plot  of  the  Ez  field.  By  visual 
inspection,  we  see  that  it  is  indeed  the  0-mode.  The  frequency  to  the  5th  decimal  can  be  found  at  the 
end  of  the  log  file  or  the  sum  file.  Record  this  value,  cf_7_o  =  34.998  ghz. 


Frequency  35.00  GHz:  Of  Ez  (V/m)  at  AREA_OUTPUT 


Another  figure  of  interest  is  the  Range  plot  of  the  Ez  field  on  Axis.  By  visual  inspection  we  see  that 
this  has  a  nice  clean  profile. 
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The  next  figure  is  the  integral  of  the  |Ez|  along  the  axis. 
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Frequency  35. 01  GHns :  Es  at  LINE  OUTPUTS 


Obtain  the  final  integrated  peak  of  Vaxis  =  9.706e6.  From  this  number  we  can  obtain  an  estimate  of 
R/Q,  using  the  following  expression. 

R/Q  =  |  V^s  |2/  (47tf  WE  )  =  1 9.706e6  |2/  (4tc  35.GHz  lJoule )  =  214. 


Note,  that  the  eigenmode  solution  is  normalized  to  We=  1  Joule. 


After  the  simulation  is  complete,  examine  the  file  “f3506.sum”.  Provided  the  eigenmode  solver 
converged  there  will  a  section  that  looks  like  the  following. 

i  ####################################################### 

I  ########  Eigenmode  Report  ############################# 

1  ####################################################### 

1  ### 

1  ###  Searching  for  a  frequency  in  the  range: 

!  ###  ...  3.37E+10  (3.50E+10)  3.63E+10 

j  ###  »>>  THE  EIGENMODE  FREQUENCY  IS:  3.50054E+10  Hz 
1  ###  ...  with  error:  0.0001  % 

1  ### 

!  ####################################################### 

1  ####################################################### 

1  ####################################################### 
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Examine  this  section  and  record  the  resultant  frequency.  If  necessary,  make  adjustments  to  the  tuners 
for  the  output  cavity.  Repeat  until  satisfied  with  the  frequency. 


Example  Step  7:  Tuning  the  O  of  the  Output  Cavity. 

In  Step  7,  we  use  a  current  driver  to  introduce  the  RF  into  the  output  cavity  at  the  frequency  of  the  0- 
mode,  much  as  we  did  for  the  single  gap  cavities,  (1  through  5).  However,  because  there  are  two 
modes  close  to  this  frequency,  we  are  not  able  to  excite  a  pure  mode.  None  the  less,  we  can  use 
measurements  similar  to  those  used  in  the  single  cavity  tuning  Steps.  .  The  section  of  the  file  for  this 
Step  looks  like  the  following: 


!****************************************************************************** 

!==>  Step7,  Output  Cavity,  Tune  Q. 

!  Cold  test  and  tune  Output  cavity. 

!  This  output  cavity  is  set  for  0-mode,  not  the  pi -mode. 

!  Measure/Tune  Q  and  R/Q 

i*********************************************** ********************* ********** 

IF  (ISTEP.eq.7)  THEN  ; 

I_Run_Type  =  ColdJTest  ; 

ColdJTest_Cavity  =  OutputCAvity  ; 

I_Cold_Type  =  2  ; 

ENDIF  ; 

PARAMETER  QT6  =  1/ (1/QE6  +I/Q06)  ;  l  Effective  Q  of  output  cavity. 

Parameter  Sigma_' Output  Cavity*  =0.0  ; 

Parameter  Sigma_'MaxCavities 1  =  0.5940E-01  ; 

!  Record  Results. 

Parameter  Q_7  =  298  ; 

Parameter  CROQ_7_RingDOWN  =  211  ; 

Parameter  CROQ_7_Eigen  =  214  ; 

The  first  figure  of  interest  is  the  measurement  for  Q.  It  will  resemble  the  following  figure.  Notice 
the  ripple  in  the  measurement.  Our  best  estimate  for  Q  is  the  center  value  between  the  extremes  (218 
minimum  to  372  maximum),  thus  we  get  Q  =  295.  There  are  additional  measurements  that  give 
similar  figures.  Verify  that  all  indicate  the  same  central  value  for  Q.  If  this  value  is  to  far  from  the 
desired  Qeffective,  than  adjust  the  conductance,  sigma  value. 
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The  next  figure  of  interest  is  the  frequency  versus  time  analysis.  Note  that  the  frequency  wobbles 
from  34.977  to  35.018  GHz,  with  center  frequency  at  34.997  GHz.  Verify  that  this  is  consistent  with 
the  eigenmode  frequency. 


Signal  analysis  Fare cpae rxc y  for  E  .  DL  sat  E _ GAP _ € _ 6 


TIME  (SEC) 


The  next  figure  of  interest  gives  us  an  estimate  of  the  R/Qaxis-  Obtain  the  value  (in  this  Step  we  have 
a  value  of  21 1.  Again  very  close  to  the  value  of  214  obtained  from  Step  6. 
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Example  Step  8;  Estimate  Mfactors  from  the  Input  Cavity. 


In  Step  8,  we  solve  cavity  1  with  eigenmode.  The  eigenmode  frequency  should  be  nearly  identical  to 
that  obtained  with  excitation  cold  test  of  Stepl .  What  we  want  from  this  simulation  is  a  measure  of 
the  Mfactor.  The  definition  we  use  is  that 

M=  M„/Md 


Where 


Mn(r)=  f  ^Z(z,r)  COS  [27lfi(z-Zcentergap)/Vbeam]  dz 

Md(r)=  I  |Ez(z,r)|  dz 


s  tin)  CE-3) 

and  the  following  figure  gives  the  profile  of  Ez(z,i=0)  cos  [2%i\ (z-Zcentergap)/ Vbeam] • 


Frequency  34 . 95  GHz  r  Ea  at?  AX  I  $ 


Using  the  integrate  option  in  MAGIC  we  obtain  the  running  integrals  of  these  profiles. 

The  maximum  value  of  Md  (r=0)  =  jEz(z,r)|  dz  =  7.485E6  is  obtained  from  the  following  figure. 


Frequency  34 . SC  GHz :  Ez  at  AXIS_CAV_1$ 


The  maximum  value  of  Mn  (r=0)  =  |  fEz(z,r)  cos  [27ifi(z-zCentergap)/vbeani]  dz  =  4.945E6  is  obtained 
from  the  following  figure. 


Frequency  34  .  SC  GHz:  Esc  at  AXIS_CAV_1$ 


e  (in)  (E— 3) 

Taking  the  ration  we  obtain  the  on-axis  M  factor  =  0.662.  Record  this  value. 

We  also  obtain  the  same  type  of  figures  at  the  beam  radius.  7.435e6 


Frequency  34.98  GHi:  Ei  at  BEAM_CAV_1$ 


Integral  value  is  5.537e6. 


Frequency  34-92  GHz :  E*  at  BEAM__CAV__1  $ 
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Taking  the  ration  we  obtain  the  on-axis  M  factor  =  0.745.  Record  this  value. 


Example  Step  9:  Estimate  Obeam  from  Short  hot  test  with  the  Input  Cavity  only. 

It  is  important  to  find  a  good  value  for  Qbeam.  This  allows  us  to  properly  set  the  desired  input  power. 

| ****************************************************************************** 
l==>  Step9.  Short  Hot  test.  Obtain  beam  loading  in  first  cavity. 

j****************************************************************************** 


IF  (ISTEP.eq.9)  THEN  ; 

I_Run__Type  =  0; 

Freq  =  CF_1  ; 
t_sim  =  5 . Onanosec ; 

ExitCavity  =  1  ; 

ENDIF  ; 

I  This  number  determined  from  hot  test.  Insert  proper  value  here. 

Parameter  Qbl  =  1206  ;  !  Beam  loading  Q. 

PARAMETER  QTbl  =  l/(l/QEl  +l/Qol  +  1/Qbl)  ; 

We  measure  the  cavityl  gap  Voltage  and  the  power  loss  in  the  empirical  load  in  cavityl.  We 
measure  Qtotal  by  evaluating  the  following  expression. 

Qtotal  =  <V2>  /  Pgap  /  (R/Q) 

We  measure  <V2>  and  Pgap  using  Magic.  We  know  the  R/Q  from  the  first  cavity. 

Once  we  have  Qtotal,  we  can  back  out  Qbeam. 

1/Qtotal  =  1/Qbeam  +  1/Qo  +  1/Qe. 

The  first  plot  of  interest  is  the  measurement  of  Pgap. 


TIME  (SEC) 
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The  next  plot  of  interest  in  this  section  is  the  measurement  of  <V2>. 


E.DL  at  L _ GAP _ 1 _ 2. 


The  next  plot  gives  us  the  estimate  of  Qtotal. 


The  next  plot  gives  us  the  estimate  of  Qbeam. 


TRANSFORM:  Q3EAM1 


The  input  power  is  given  by  the  following  figure. 


Example  Step  10:  Run  Frequency  and  Power  Sets. 

We  use  the  following  expression  to  estimate  the  necessary  current  for  the  desired  input  power.  We 
first  get  the  drive  power  required  from  the  following  expression. 


P Drive  PfDesired  /(I  —  Qtbl/Qel) 


Since  we  the  driver  is  actually  a  current  source  we  need  to  relate  the  applied  current  to  the  applied 
power.  We  note  that  the  power  can  be  written  at  the  applied  current  times  the  gap  voltage.  We 
estimate  the  gap  voltage  from  the  following  expression. 

Vgap  =  Igap  Q  (R/Q)  l/[l+i2Q(f-f,)/f]  Mfactor 

Define  a  phase  factor,  <f>,  with  the  following  expression: 

i/[i+i2Q(f-fi)/f]  | 


13 


Then  we  finally  use  the  following  expression  to  set  the  applied  current. 


Idrive  (2  Pdrive  /  (Qtbl  (R/Q)l  $  Mfact0r  ) 

The  only  user  inputs  for  Step  10  are  the  frequency  and  the  desired  input  power. 

For  this  design  the  range  of  frequencies  recommended  is  34.8  GHz  to  35.2  GHz,  and  the  recommend 
input  power  is  0.1  milliwatt  to  100  milliwatt. 

J  *********************************************************************** 

!==>  CaselO.  Hot  tests.  Select  frequency  and  input  power. 

i*********************************************************************** 

IF  (istep.eq.10)  THEN  ; 

I_Run_Type  =  0 ; 

t_sim  =  50.0nanosec; 

Header  Remarks  "Hot  Test  "  ; 

!  Select  Drive  Frequency  and  Input  Power  for  Hot  tests. 
parameter_ghz  Freq  =  35.000  GHz; 

Parameter  Desired_PowerIn  =  l_milliwatts; 

!  These  are  the  nominal  operating  parameters. 

ParameterJfCiloVolt  V_beam  =  14.0kV  ; 

Parameter_Amp  A_beam  =  0 . 5A  ; 

Parameter_Gauss  Bz  -  0.80tesla; 

ENDIF  ; 

The  following  figures  illustrate  the  type  of  output  available  for  each  hot  test  simulation. 

The  simulation  below  was  at  35  Ghz  with  an  input  power  of  0.1  milliwatt. 

Input  power:  0.1  mWatt 


TRANSFORM:  POWEREXT 


Volts  DB  (  Volts/GHz) 


FFT  of  Voltage  in  output  cavity.  1.26kV 


Magnitude  Of  FFT  Of  E.DL  at  L_GAP_6_7 


E.DL  at  L  GAP  1  1 


0  10  20  30  40  SO 

TIME  (SEC)  (E- 9) 


Voltage  in  cavity  2;  10.44 

E.DL  at  L  GAP  2  2 


o 

N 


O 


a 

N 

I 


I 

20 

TIME  (SEC) 


50 

(E-9) 


0 


10 


-r 

30 


-h 

40 


APPENDIX:  Input  file  listing 
F3506.m2d 


I  ********************************************************************** 

!  F3607 . m2d  -  Base  File  Name  for  35  GHz,  7  Cavity  klystron. 

I  ********************************************************************** 

!  Cavity  1:  Input  Cavity. 

!  Cavity  2-5:  Intermediate  Cavities. 

!  Cavity  6&7:  Dual  Output  Cavity 

!  Cold  test  runs  of  individual  cavities.  Excite  using  antenna  driver  at 
!  the  target  cavity  frequency.  Adjust  by  adding/ removing  conductor  cells 
!  from  the  the  tuning  strips.  Resolution  is  approximately  6.5Mhz/cell. 

!  Adding /removing  cells  shifts  frequency  up/down  as  the  inductance 
!  changes. 

Comment  "  Specify  output  cavity.  " ; 

Integer  OutputCavity,  ExitCavity,  ColdTestCavity  ; 

Output Cavity  =  6  ; 

MaxCavities  =  Output Cavity+1  ; 

ExitCavity  =  MaxCavities  ; 

ColdTestCavity  =  1; 

Comment  "  The  desired  Cavity  Design  frequencies."; 


PARAMETERJ3HZ 

f  tube 

= 

35.0GHz; 

l  operating  mode  of 

the  output 

PARAMETER_GHZ 

fl 

= 

34.98  GHz 

PARAMETERJ3HZ 

f  2 

= 

34.90  GHz 

!  " 

PARAMETER__GHZ 

f  3 

= 

35.09  GHz 

!  11 

PARAMETER_GHZ 

f  4 

= 

35.09  GHz 

i  « 

PARAMETER_GHZ 

f  5 

= 

35.17  GHz 

1  " 

PARAMETER_GHZ 

f  6 

= 

35.00  GHz 

I  " 

PARAMETER_GHZ 

fOut_0 

= 

34.999  GHz; 

!  Use 

Eigenmode  Dbl 

prec . 

PARAMETERJ3HZ 

fOut_PI 

= 

34. 4564GHz; 

!  Use 

Eigenmode  Dbl 

prec . 

Comment  "  Nominal  operating 

parameters . 

H  . 

t 

parameterjjhz 

Freq 

35.000 

GHz; 

Parameter  Desired  Powerln  =  1  milliwatts; 


Parameter_KiloVolt  V_beam 
Parameter_Amp  A__beam 

Parameter_Gauss  Bz 

Parametert  t  sim 


14  .  OkV  ; 

0.5A  ; 
0.80tesla; 

40 . Onanosec; 


Comment  "  The  desired  Cavity  Qs  and  R/Q."; 


Parameter 

QE1 

= 

350 

Parameter 

Qol 

= 

2000 

Parameter 

Rover Q1 

= 

128 

Parameter 

QE2 

= 

1E9 

Parameter 

Q02 

= 

2000 

Parameter 

RoverQ2 

= 

128 

Parameter 

QE3 

= 

1E9 

Parameter 

Q03 

= 

2000 

Parameter 

RoverQ3 

= 

128 

Parameter 

QE4 

= 

1E9 

Parameter 

Q04 

= 

2000 

Parameter 

RoverQ4 

= 

128 

Parameter 

QE5 

= 

1E9 

Parameter 

005 

= 

2000 

Parameter 

RoverQS 

= 

128 

Parameter 

QE6 

= 

350 

Parameter 

Q06 

2000 

Parameter 

RoverQ6 

= 

128 

Comment  "  Design  Gap  to  Gap  spacing. 


P  ARAME  TE  R_MM 
PARAMETER_MM 
PARAMETER_MM 
P ARAME  TE  R_MM 
PARAMETER_MM 
PARAMETER__MM 
PARAMETERJ4M 
PARAMETER  MM 


dZ_gap2gap_0_l 

dZ_gap2gap_l_2 

dZ_gap2gap_2_3 

dZ_gap2gap__3_4 

dzL9aP29aP_4_5 
dZ_gap2gap__5_6 
dZ_gap2  gap_6__7 
dZ__gap2gap_7_8 


5.0  mm; 
5.0  mm; 
5.0  mm; 
5.0  mm; 
5.0  mm; 
4.7  mm ; 
2.0  mm; 
5.0  mm; 


!  Parameter  definitions.  Do  NOT  change  these. 


On  «  1  ; 

Off  a  0; 

Hot_Test  =  0; 

Cold_Test  =  1 ; 

EigenMode_Test  =  2; 

I_Cav_Cndctnc  =  on  ;!  Use  cavity  loss  in  all  cavities  if  present. 
I_Run_Type  =0;  10=  Hot  test 

!  1  =  Antenna  Driven  cold  test 
!  2  =  Eigenmode  cold  test 

I_Cold_Type  =2  ;  !  1  =  lbyl  (applies  to  all  cavities  and  cells) 

~  !  2  =  2by2  (applies  only  to  output  structure) 

!  3  =  2by2  pulse 

l  Cold  test  cavity,  or  first  cell  if  I_Cold__Type  =  2. 
ColdJTest_Cavity  =  OutputCAvity  ; 

!  Header  Definition. 


HEADER  AUTHOR  "Liqun  Song"  ; 

HEADER  DEVICE  "35GHz  Klystron"; 

HEADER  ORGANIZATION  "SLAC"; 

GRAPHICS  PAUSEon; 
graphics  pause  ; 

!  Run  Controls  and  Global  Parameters  ******************* 
!  Select  simulation  cases  here. 


IStepal 

IStep=2 

IStep=3 

IStep=4 

IStep=5 

IStep=6 

IStep=7 

IStep=8 

IStep=9 

IStep=10 


!  Magic2d__Sng:  Cold  test,  Tune  cavity  1. 

!  Magic2d_J3ng:  Cold  test,  Tune  cavity  2. 

!  Magic2d_Sng:  Cold  test,  Tune  cavity  3. 

!  Magic2d__Sng:  Cold  test,  Tune  cavity  4. 

I  Magic2dJ3ng:  Cold  test.  Tune  cavity  5. 

!  Magic2d_Dbl :  Eigenmode  test,  Tune  cavity  6  Frequency. 

!  Magic2d_Sng:  Cold  test,  Tune  cavity  6,  Tune  Q. 

!  Magic2d__Dbl :  Eigenmode  cavity  1.  Measure  M  factors. 

1  Magic2d_Sng:  Hot  test  of  cavity  1,  obtain  Beam  loading  Q. 
!  Magic2d_Sng:  Run  Hot  tests  at  desired  frequency  and 
!  input  power. 


t  ***  select  case  here. 
ISTEP  =  1  ; 


*1' 


* 


I*********************************************************************** 

i  Cavity  Tuning  Procedure,  for  Input  and  Intermediate  Cavities. 

!  Procedure,  run  Case  and  obtain  frequency. 

!  If  frequency  lo,  increase  values  of  tune_left  and  tune_right . 

!  If  frequency  hi,  decrease  values  of  tune_left  and  tune_right . 

!  IF  frequency  good.  Check  Qeffective  and  R/Q. 

!  Remember  that  you  are  measuring  and  tuning  to  Q  effective 
1  IF  Q  (measured)  is  greater/less  than  Q  (desired)  , 

!  then  alter  value  of  sigma. 

!  Adjust  sigma  as  follows: 

!  Sigma  (desired)  =  Sigma  (used)  *  Q  (measured)  /  Q  (desired)  . 

i*********************************************************************** 

!==>  Casel,  Cavityl,  Can  Use  either  Magic2d_Sng  or  Magic2d_Dbl 
!  Cold  test  and  tune  cavity  1. 

j *********************************************************************** 


IF  (IStep.eq.l)  THEN  ; 

I_Run_Type  =1;  !  0  =  Hot  test 

Cold__Test__Cavity  =  1  ; 

ExitCavity  =  1  ; 

Header  Remarks  "Tuning  Input  Cavity"  ; 

ENDIF  ; 


!  Cavity  1,  tuning  parameters. 

Parameter__GHZ  Desired_fl  =  fl;  ! 
Tune_Left_l=  13  ;  I 
Tune__Right_l=  12;  ! 


PARAMETER  QT1  -  l/(l/QEl  +l/Qol  ) 

Parameter  Sigma_l  =  0.32114E-01; 

!  Record  results  here. 
ParameterJSHZ  CFJL  =  34.982  GHZ  ; 

Parameter  CQ_1  =  298  ; 

Parameter  CR0Q_1  =  127 


Desired  Frequency  for  Cavity  1 
Use  tuners  to  adjust  the  frequency, 
Change  by  about  +/-1,  to  give 
about  +/1  6.5  MHz  frequency  change. 

;  !  Desired  Q  effective  (of  Qo  and  Qe) 

!  Change  this  value  to  adjust  the  Q. 


!  Frequency 

I  Q  effective  (of  Qo  and  Qe) 

!  R  over  Q  as  measured  on  Axis. 


!==>  Case2,  Cavity2,  Can  Use  either  Magic2d_Sng  or  Magic2dJDbl 
I  Cold  test  and  tune  cavity  2 . 


IF  (IStep.eq.2)  THEN  ; 

I_Run_Type  =  1  ; 

Cold_Test_Cavity  =  2  ; 

ExitCavity  =  2  ; 

Header  Remarks  "Tuning  Intermediate  Cavity  2"  ; 
ENDIF  ; 


Parameter_GHZ  Desired_f2  =  f2; 
Tune_Lef t_2=  6  ; 

Tune_Right_2=  6  ; 


!  Desired  Frequency  for  Cavity  2 
I  Use  tuners  to  adjust  the  frequency, 
!  Change  by  about  +/-1,  to  give 
!  about  +/1  6.5  MHz  frequency  change. 


PARAMETER  QT2  =  1/(1/QE2  +l/Qo2) 
Parameter  Sigma_2  =  0.47441E-02; 


1  Desired  Q  effective  (of  Qo  and  Qe) 
!  Change  this  value  to  adjust  the  Q 


!  Record  results  here. 
Parameter_GHZ  CF_2  =  34.898  GHZ  ; 
Parameter  CQ_2  =  2000; 
Parameter  CR0Q_2  =  127  ; 


ai 


J  ******************************************************************************** 

1 =->  Cases,  Cavity3,  Can  Use  either  Magic2d_Sng  or  Magic2d_Dbl 
!  Cold  test  and  tune  cavity  3, 

j ******************************************************************************** 

IF  (IStep.eq. 3)  THEN  ; 

I_Run__Type  =  1 ; 

Cold__Test_Cavity  =  3  ; 

ExitCavity  =  3  ; 

Header  Remarks  "Tuning  Intermediate  Cavity  3"  ; 

ENDIF  ; 

Parameter_GHZ  Desired_f3  =  f3;  !  Desired  Frequency  for  Cavity 

Tune_Left_3=  21  ;  !  Use  tuners  to  adjust  the  frequency, 

Tune_Right_3=  21;  !  Change  by  about  +/-1,  to  give 

!  about  +/1  6.5  MHz  frequency  change. 

PARAMETER  QT3  =  l/(l/QE3  +l/Qo3)  ;  !  Desired  Q  effective  (of  Qo  and  Qe) Parameter  Sigma_3 

=  0 . 48437E-02 ;  !  Change  this  value  to  adjust  the  Q. 

!  Record  results  here. 

Parameter_GHZ  CF_3  =  35.091  GHZ  ; 

Parameter  CQ_3  =  2000; 

Parameter  CR0Q_3  =  127  ; 

I******************************************************************************** 

! ==>  Case4,  Cavity4,  Can  Use  either  Magic2d__Sng  or  Magic2d_Dbl 
!  Cold  test  and  tune  cavity  4 . 

1******************************************************************************** 

IF  (IStep.eq. 4)  THEN  ; 

I_Run_Type  =  1 ; 

Cold_Test_Cavity  =  4  ; 

ExitCavity  =  4  ; 

Header  Remarks  "Tuning  Intermediate  Cavity  4"  ; 

ENDIF  ; 

Parameter_GHZ  Desired_f4  =  f4;  !  Desired  Frequency  for  Cavity 

Tune__Left_4  =  21  ;  !  Use  tuners  to  adjust  the  frequency, 

Tune_Right_4=  21;  !  Change  by  about  +/-1,  to  give 

!  about  +/1  6.5  MHz  frequency  change. 

PARAMETER  QT4  =  l/(l/QE4  +l/Qo4)  ;  !  Desired  Q  effective  (of  Qo  and  Qe) 

Parameter  Sigma__4  =0.48437E-02;  !  Change  this  value  to  adjust  the  Q. 

I  Record  results  here. 

Parameter_GHZ  CF_4  =  35.091  GHZ  ; 

Parameter  CQ__4  =  2000; 

Parameter  CR0Q_4  =127; 

j ******************************************************************************** 

!==>  Case5,  CavityS,  Can  Use  either  Magic2d_Sng  or  Magic2d__Dbl 
!  Cold  test  and  tune  cavity  5. 

j ******************************************************************************** 

IF  (IStep.eq. 5)  THEN  ; 

I_Run_Type  =  1; 

Cold_Test_Cavity  =  5  ; 

ExitCavity  =  5  ; 

Header  Remarks  "Tuning  Intermediate  Cavity  5"  ; 

ENDIF  ; 

Parameter_GHZ  Desired__f5  =  f5;  !  Desired  Frequency  for  Cavity 


a,  o  r 


Tune_Le  f  t_5 =  2  ; 
Tune_Right_5=  3 ; 


!  Use  tuners  to  adjust  the  frequency, 
!  Change  by  about  +/-1,  to  give 
!  about  +/1  6.5  MHz  frequency  change. 


PARAMETER  QT5  =  1/ (1/QE5  +l/Qo5)  ;  !  Desired  Q  effective  (of  Qo  and  Qe) 

Parameter  Sigma_5  =  0.48861E-02;  !  Change  this  value  to  adjust  the  Q. 


!  Record  Results. 

Parameter_GHZ  CF_5  =  35.168  GHZ  ; 
Parameter  CQ_5  =  1996; 

Parameter  CR0Q_5  =  127  ; 


t******************************************************************************** 

! ==>  Case6,  Output  Cavity,  Use  Magic2d_Dbl 
!  Cold  test  and  tune  Output  cavity. 

!  This  output  cavity  is  set  for  0-mode,  not  the  pi -mode. 

!  Output  Cavities,  USe  double  precision  Eigenmode  to  tune  these  cavities. 

!  Get  R/Q  FROM  EIGENMODE,  Integrate  |Ez|  along  axis. 

I  Then  R/Q  =  ( |  E  |  .dl)  A2/4pi/f /ljoule__f  rom  eigen  soln. 
!******************************************************************************** 


IF  (IStep.eq.6)  THEN  ; 

I_Run_Type  =  2  ; 

Cold__Test_Cavity  =  OutputCAvity  ; 

I_Cold_Type  =  2  ; 

Header  Remarks  "Frequency  Tuning  Output  Cavity  "  ; 

ENDIF  ; 

Parameter_GHZ  Eigenmode_f0ut_0  =  Fout_0  ; 

Tune_Le  f  t_6 =  8 ; 

Tune_Right_6=  8; 

Tune_Le  f  t_7 =  8 ; 

Tune_Right_7=  7; 

i  Record  Results. 

Parameter_GHZ  CF_7__0  =  35.0054  GHZ;  !  adjusted  after  eigenmode. 


t******************************************************************************** 

I==>  Case7,  Output  Cavity. 

!  Cold  test  and  tune  Output  cavity. 

!  This  output  cavity  is  set  for  0-mode,  not  the  pi-mode. 

!  Measure/Tune  Q  and  R/Q 

I******************************************************************************** 


IF  (IStep.eq.7)  THEN  ; 

I_Run_Type  =  ColdJTest  ; 

Cold_Test_Cavity  =  OutputCAvity  ; 

I_Cold_Type  =  2  ; 

Header  Remarks  "Tuning  Q/Load  of  Output  Cavity  "  ; 

ENDIF  ; 

PARAMETER  QT6  =  1/(1/QE6  +I/Q06)  ;  !  Effective  Q  of  output  cavity. 

Parameter  Sigma_' Output Cavity'  =  0.0  ; 

Parameter  Sigma_' MaxCavities '  =  0.5940E-01  ; 

!  Record  Results. 

Parameter  Q_7  =  298  ; 

Parameter  CROQ_7_R i ngDO WN  =  211  ; 

Parameter  CROQ_7_Eigen  =214  ; 


I******************************************************************************** 

!  ==>  Case8,  Cavityl,  Can  Use  either  Magic2d_Sng  or  Magic2d_Dbl 
I  Eigenmode  test  and  tune  cavity  1. 


************************************************************************ 
IF  (IStep.eq. 8)  THEN  ; 

I_Run_Type  =  2  ; 

Coid_Test_Cavity  =  1  ; 

Freq  =  CF_1  ; 

ExitCavity  =  1  ; 

Header  Remarks  "Cavityl,  obtain  Beam  Mf actors  "  ; 

ENDIF  ; 

Parameter  KiloVolt  V  beam  =  14.0kV  ; 


!  Record  Results. 

Parameter  Axial_M_f actor  =  0.662  ;  !  these  are  obtained  from  range  integrals, 

!  AND  eigenmode  OF  CAVITY1. 

Parameter  BeamR_M_f actor  =  0.745  ; 

i *********************************************************************** 

!==>  Case9 .  Short  Hot  test.  Obtain  beam  loading  in  first  cavity. 

I*********************************************************************** 

IF  (IStep.eq. 9)  THEN  ; 

I_Run_Type  =  0  ; 

Freq  =  CF_1  ; 
t_sim  =  12.0nanosec; 

ExitCavity  =  1  ; 

Header  Remarks  "Cavityl,  obtain  Beam  Loading  Q  "  ; 


!  Select  Drive  Frequency  and  Input  Power  for  Hot  tests. 


parameter_ghz  Freq 

Parameter  Desired_PowerIn 
Parameter_KiloVolt  V_beam 
Parameter_Amp  A_beam 

Parameter_Gauss  Bz 

ENDIF  ; 


35.000  GHz; 
l_milliwatts ; 
14 . OkV  ; 

0.5A  ; 
0.80tesla; 


!  This  number  determined  from  hot  test. 

!  Record  proper  value  here. 

Parameter  Qbl  =  1206  ;  !  Beam  loading  Q. 

PARAMETER  QTbl  =  1/(1/QE1  +l/Qol  +  1/Qbl)  ; 

i************************************************************** ********* 

!==>  CaselO.  Hot  tests.  Select  frequency  and  input  power. 

j  *********************************************************************** 


IF  (IStep.eq. 10)  THEN  ; 
I__Run_Type  =  0  ; 

t_sim  =  BO.Onanosec; 
Header  Remarks  "Hot  Test  " 


I  Select  Drive  Frequency  and  Input  Power  for  Hot  tests. 


parameter_ghz  Freq 

Parameter  Desired_PowerIn 
Parameter_KiloVolt  V_beam 
Parameter_Amp  A_beam 

Parameter_Gauss  Bz 

ENDIF  ; 


35.000  GHz; 
l__milliwatts; 
14. OkV  ; 

0.5A  ; 

0 . 80tesla; 


Start  internal  template. 


TEMPLATE  KLYSTRON  ; 
Stop  ; 


